Andreev reflection through Fano resonances in molecular wires 
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We study Andreev reflection in a normal conductor-molecule-superconductor junction using a 
first principles approach. In particular, we focus on a family of molecules consisting of a molecular 
backbone and a weakly coupled side group. We show that the presence of the side group can lead 
to a Fano resonance in the Andreev reflection. We use a simple theoretical model to explain the 
results of the numerical calculations and to make predictions about the possible sub-gap resonance 
structures in the Andreev reflection coefficient. 

PACS numbers: 73.63.-b,74.45.+c 



Fano resonances^ are a universal interference phe- 
nomenon which can affect coherent electrical trans- 
port through nanostructures in many different sys- 
tems. Examples of Fano lineshape in mesoscopic sys- 
tems include scanning tunnelling microscope measure- 
ments on a single magnetic atom absorbed on a gold 
surfaceSii^, single-electron transistors fabricated into a 
gated two-dimensional electron gas^, quantum dots 
embedded into an Aharonov-Bohm ring^i^, multiwall 
carbon nanotubesZi&iS and recently single-wall carbon 
nanotubes^ and double-wall nanotubes^i. Fano reso- 
nances (FRs) also appear in the the conductance of quasi 
one-dimensional quantum wires with donor impuritiesi^ 
and in the case of quantum wires with a side coupled 
quantum doti^. In molecular electronics, due to the re- 
alistic treatment of the metal electrodes, FRs have been 
found in the transmission of dithiol benzeneii>. More 
generally, theoretical calculations predict that Fano- 
lineshape should appear in the transmission through 
molecular wires with attached side groups^^ or as a con- 
sequence of quantum interference between surface states 
of the measuring electrodes and the molecular orbitala^^. 

If one of the measuring probes is superconducting, 
the conductance for energies E smaller than the su- 
perconducting pair potential A depends on the An- 
dreev reflection probability Ra{E). The Andreev re- 
flection in various mesoscopic systems has been stud- 
ied for a long time (see e.g Refs. llTlflsl and references 
therein) but the interest has recently renewed when An- 
dreev reflection through carbon nanotubes was measured 
experinientallyi^i^. These experiments have sparked nu- 
merous theoretical studies both in the absence of the 

electron-electron interactio n^^i^^ and in the presence of 
the interactioii2^i25;26^j28^^^^^ j^^ ^^^^ ^j^ggg 

studies it was assumed that it was sufRcient to con- 
sider resonant transport through a single energy level 
and as a consequence, the Andreev reflection as a func- 
tion of energy exhibited Breit-Wigner type resonances. 
A notable exception is Ref. 'ss' where transport through 
an Aharonov-Bohm ring with an interacting quantum 
dot situated in one of its arms was considered and a 
Fano-type asymmetric resonance was found in the con- 
ductance. Very recently, Tanaka et al.'^^ studied An- 
dreev transport through side-coupled interacting quan- 
tum dots focusing on the interplay of Andreev scattering 



and Kondo effect. 

It was demonstrated in Ref. [l^ that Fano resonances 
are a generic feature of molecular wires with attached 
side groups. It was also shown that for a certain type of 
molecular wires a FR can appear in the normal conduc- 
tance Gn{E) very close to the Fermi energy Ep. In a 
normal metal-molecule-superconductor (N-Mol-S) junc- 
tion therefore these FRs would also affect the sub-gap 
transport. The aim of this paper is to study Andreev 
reflection through molecular wires when the normal con- 
ductance exhibits FRs close to the Fermi energy. Per- 
forming ab initio simulations of molecular wires in N- 
Mol-S junctions we show how FRs influence the sub-gap 
transport. We elucidate the results of the numerical cal- 
culations using a simple analytic model. We also predict 
that for finite energies the differential conductance can 
reach the unitary limit if there is a strong asymmetry in 
the coupling to the leads. 



Au 



Au 



FIG. 1: Possible experimental setup. The molecule is con- 
tacted by gold electrodes, one of which is superconducting 
due to the proximity effect. 

A possible experimental setup to measure Andreev re- 
flection in N-Mol-S junctions is shown in Fig. [1] The 
molecule is contacted with thin gold electrodes on both 
sides. On top of one of the electrodes a second layer of 
e.g. aluminium or niobium is deposited, which at low 
enough temperature becomes superconducting. Due to 
the proximity effect this top layer induces superconduc- 
tivity in the gold electrode beneath (in our calculations 
we assume that the induced superconductivity is s type). 
We note that this setup was successfully used in Ref. [20 
to study Andreev-reflecion in normal conductor - carbon 
nanotube - superconductor (N-Cn-S) junctions. 

To study FRs in a N-Mol-S system, we choose the 
smallest molecule of a recently synthesized family of 
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molecular wires^'^'^. Since these molecules have ter- 
minal thiol groups they can easily bind to gold surfaces 
making them ideal for experiments on single-molecule 
transport properties. The central part of the molecule 
consists of a single fluorenone unit, which could be chem- 
ically modified, e.g. by replacing the oxygen with bipyri- 
dine rings, as shown schematically in Fig. [T] The differ- 
ential conductance of the system was calculated using a 
combination of the DFT code SIESTA^^^ and a Green's 
function scattering approach explained in Refs. [39ll40l . 
Initially the isolated molecule is relaxed to find the op- 
timum geometry, then the molecule is extended to in- 
clude surface layers of the gold leads. In this way, charge 
transfer at the gold-molecule interface is included self- 
consistently. The number Ng of gold layers is increased 
until computed transport properties between the (nor- 
mal conducting) gold leads no longer changed with in- 
creasing Ng. Typically, this extended molecule contained 
A'g = 3 to 4 gold layers on each side, and the layers con- 
sisted of 9 atoms on the (111) plane. The leads, which 
were assumed to be periodic in the transport direction, 
also consisted of gold layers containing 9 atoms on the 
(111) plain. Using a double-C basis plus polarization 
orbitals, TrouUier-Martins pseudopotentials^ and the 
Ceperley-Alder LDA method to describe the exchange 
correlationiS^ effective tight-binding Hamiltonians Hm, 
Hl of the extended molecule and of the leads, respec- 
tively, were obtained. To investigate the generic physics 
of this system, we employ the simplest possible approxi- 
mation for the order parameter, namely that it changes 
in a step-function-like manner at the superconducting 
lead - extended molecule interface. Therefore the super- 
conducting lead was modelled by introducing couplings 
of constant magnitude A between the electron and hole 
degrees of freedom in H^, while no such coupling was 
present in the extended molecule and in the normal lead. 
We focus on sub-gap transport and therefore compute the 
Andreev- reflection probability Ra{E), because at zero 
temperature for E < A the differential conductance is 
given by Gns{E) = ^Ra{E). 

The calculations of Ref. [l^ have shown that by chang- 
ing by rotational conformation of the bipyridine unit 
it is possible to change the position of the Fano reso- 
nance with respect to Ep. The definition of the angle 
of rotation 9 of the bipyridine group is the following: 
9 = Q° when the rings of the sidegroup are parallel to the 
molecule axis and it is 90° when they lie perpendicular. 

We consider the molecule whose rotational conforma- 
tion is 6* = 71.4°. Assuming first that both leads are nor- 
mal conducting (N-Mol-N junction), close to the Fermi 
energy there is a FR in the differential conductance 
Gn{E) = ^Tn{E) (where Tn{E) is the normal trans- 
mission) as it can be seen in Fig. [21 Since for conventional 
supercondcutors the typical superconducting gap values 
are 0.1 — 1.5 meV, we first consider the transport for en- 
ergies i? < A = 1.35 meV. The Fano peak in Gn{E) is 
at 6E w 3.7meV above Ep, therefore SE is bigger than 
A and the infiuence of this resonance on the sub-gap 



transport can be understood by considering the zero bias 
conductance Gns{0)- Indeed, as Fig. [3] shows, the An- 
dreev reflection is almost constant apart from the region 
E A where a sharp peak can be observed which is due 
to the singularity in the density of states of the super- 
conductor at this energy. One can see that off-resonance 




5 10 15 20 25 30 35 40 45 



E-Ep [mcV] 

FIG. 2: Normal differential conductance Gjv(-E) (in units of 
2e^/h, dashed line) and sub-gap conductance Gns{E) (in 
units of 4e^//i, solid line) in logarithmic scale as a function of 
energy. We used A = 41 meV. 

Gns{0) is smaller than Gjv(O). 

More generally however, if there is a narrow F < A 
resonance at some \5E\ < A above or below Ep, the en- 
ergy dependence of the Andreev reflection becomes im- 
portant. (In case e.g. of carbon nanotubes, which can 
be gated, this scenario should be easily attainable, as in 
Ref 10 where the width of the Fano peak was ~ 0.2 mcV.) 
To illustrate this case, we performed computations using 
the same molecule but much bigger A. The results of 
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FIG. 3: TnIE) (dashed line) and RAiE) (sohd line) as a 
function of energy. We used A = 1.35 meV. The sharp peak 
in Ra{E) is at iJ ~ a (see main text). 

the computations are shown in Fig. [51 As one can see, a 
Fano resonance now appears both Gat in Gns- However, 
a closer inspection reveals that the width of the Fano 
peak in Gns at « 3.7 meV is roughly half of the width 
of the corresponding peak in Gat. 

To explain the results of the numerical calculations we 
consider a simple model, introduced in Ref. 15, which was 
shown to capture the essential features of the transport 
between normal conducting leads. Close to a resonance. 
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it is sufRcient to consider a single backbone state 
with resonant energy ii and a state I/2) of energy £2 
which is associated with a side group of the molecule 
(ei and £2 are measured relative to the Ep). The weak 
coupling between the backbone of the molecule and the 
side group is described by a matrix H12. We denote by 
tc = {f2\Hi2\fi) the couphng between the two states, 
whereas the coupling of the backbone state to the normal 
(superconductor) lead is described by matrices Wjsr (Ws)- 
A brief derivation of the Andreev reflection probability 
Ra{E) for this system is given in Appendix 1X1 here we 
only summarize the main results. 
The linear conductance is given by 



Gns{0) 



4e2 



4p2 p2 =i 



h [ie+s-r + irj 



(1) 



Here Tl, (Th) is the normal state tunnelling rate to the 
left (right) lead at Ep and e± = £ ± V 6e^ + where £ = 
(£1 +£2)72, Se — {si — £2)72. The maximal conductance 
is attained at Tl = ^r, £± = when it is twice as large 
as the normal conductance. Note that the conductance 
maximum is not attained when £1 is aligned with Ep 
as one might expect. The hybridization between £1 and 
£2 due to the coupling tc leads to a different resonance 
condition for this system. Off-resonance, i.e. when £-1- ^ 
, Gns falls off more rapidly as a function of £± than Gn 
[see Eq. (1) in Ref . ITsj . Therefore the G'Ar5(0) is usually 
smaller than Gjv(O). Moreover, Gns is zero if £2 = 0, 
i.e. when the energy of the side coupled state equals Ep. 
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FIG. 4: Gn (in units of 2e^ /h, dashed) and Gns (in units 
of Ae^ /h, solid) in logarithmic scale as a function of energy. 
We used \/T%r%/\af'\ = 4 in the case of a), b), c) and 
= 0.25, in the case of d). 

For finite energies < A the most important features 
of the differential conductance of our model are the fol- 
lowing. If the coupling to the normal lead is stronger than 

> 



to the superconducting one, i.e. when 
where (^n) tunnelling rates for electrons (holes) 



pe p/i 
^ AT 



from the normal lead and cr|'* is an off-diagonal element 
of the self-energy matrix S5 (see Appendix \X\ for the pre- 
cise definitions, as well as for the definitions of e'^, £l, 
to be introduced below), in good approximation 



4e2 

Gns{E)^—A{E)Tc{E) 



(2) 



where the amplitude A{E) is a slowly varying function 
of the energy and 



Te{E) 



(3) 



[^E~el){E-etW + {T%riE~S2r 

Assuming a weak coupling between |/i) and I/2) i.e. that 
tc <C fe = |£+ — £f.| (which also means that £f_ ^ £2), 
for energies close to £l the probability amplitude Te{E) 
can be further approximated by 



Tc{E)^A 



a^(p- -I- 1 



(4) 



where A = r%T%/{et - e^f ., e ^ {E ~ e<t_)lT 
(£l - e+)7(£i - e2f and q = {et - e2)lV%. There- 
fore, if < £'L,£2 < A a FR will appear in the sub- 
gap transport [see FigdJ^a)]. For strong couphng such 
that ^ \et_ — £2! the Fano lineshape would become 
a symmetric dip. If however £^ < A < £'L,£2 is satis- 
fied, a Breit-Wigner resonance (BWR) of width F^ oc- 
curs [shown in Fig. HKb)], while for < £!j_,£l,£2 < A 
the Andreev reflection exhibits both a FR and a BWR 
[Fig.HI^c)]. Note, that Te{E) is very similar to the trans- 
mission amplitude Tn{E) calculated in Ref. [l^ for nor- 



N- 



pe p/i 



> 



the 



mal conducting leads. Since for 

resonance energies £f_, e'^ are usually very close to the 
resonance energies appearing in the expression of r/v(i?), 
one finds that the resonance structures of the normal 
conductance will also appear in the sub-gap transport 
if the relevant resonance energies are smaller than the 
superconducting pair potential. This explains the occur- 
rence of a Fano resonance in Gns{E) in Fig. (5] How- 
ever, since A{E) in Eq. ^ is usually much smaller than 
unity, G]ys{E) itself can also be smaller than Gm{E). 
The widths of the resonances in the Andreev-reflection 
coefficient can be significantly smaller than in the nor- 
mal transmission. This happens because coupling to the 
superconductor does not lead to the broadening of the 

resonant levels. Therefore if 



^.''l the peaks 



in the normal and in the Andreev transport have roughly 
the same width while for \Jt%X^ > \af\ the width of 

the peaks in the Andreev reflection is half of the width 
of the corresponding peaks in the normal transmission. 
This can also be observed in Fig. ^ We note that for 
E fa A where a^^ changes rapidly with energy the for- 
mula shown in Eq. ^ is not applicable because in the 
derivation of Eq. ^ we have assumed that the self energy 
CTc/* is a slowly varying function of the energy. 
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Finally, we briefly discuss the predictions of our model 
for the case when the coupling to the superconductor 
is stronger than to the normal lead, i.e. when Icrf'*! > 

y'r^r^. The conductance can no longer be approxi- 
mated by Eq. ([2]) because a^^ introduces hybridization 
between electron and hole levels. We find that in the 
most general case the conductance exhibits both a FR 
and a BWR, if the corresponding resonance energies are 
smaller than the superconducting gap. These peaks, as 
mentioned before, can be much narrower than the ones 
in the normal transmission because the superconductor 
docs not broaden them. Moreover, we find that for 



> 



r^r^ the conductance can even reach the 



unitarity limit. This could not happen in the opposite, 

I I 



case because a resonance in Te{E) is 

not accompanied by a resonance in A{E) and therefore 
the conductance is always smaller than Ae'^/h. We illus- 
trate this in Fig. [UJd) where Gns is shown along with 
Gn- One can see that Gn < 2e^//i because the cou- 
plings to the leads are asymmetric and there is a broad 
resonance at E/ A « 0.15 along with an almost symmet- 
ric, narrow dip at E/A « 0.32. In contrast, Gns has a 
narrow FR and also a BWR, the latter peak reaching the 
unitarity limit. 

In summary, we have studied the Andreev reflection 
through a class of molecules which exhibit Fano reso- 
nances in the normal conductance. Our numerical cal- 
culations based on ab initio methods indicate that Fano 
resonances may also appear in the sub-gap transport. A 
simple theoretical model that we used to understand the 
results of the numerical calculations predicts that a) if 
the coupling to the normal lead is weaker than the cou- 
pling to the superconducting one, the resonance structure 
of the normal conductance can manifest itself in the An- 
dreev reflection coefficient if the resonance energies are 
smaller than the superconducting gap and b) if the cou- 
pling to the superconductor is strong, the resonances in 
the normal conductance and in the Andreev reflection 
can be very different, both in position and in width. 
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superconductor hybrid systems^. Here we employ the 
Green's function technique presented in Ref. [1^ in which 
the Hilbert space is divided into a sub-space A contain- 
ing the external leads and a sub-space B containing the 
molecule. 

Assuming for a moment that the molecule is isolated, 
for energies close to a resonance it can be described by 
quantum states I/2) with resonant energies ei, £2- 
These states are coupled together by a hamiltonian H12 
with matrix element tc = (/i|-ffi2|/2)- The effect of cou- 
pling of the molecule to the normal conducting (super- 
conducting) lead via a coupling matrix Wn (Ws) is rep- 
resented by the energy dependent self-energy matrices 
Sjv = cTjv - iTpf (Ss = (Ts - iTs) where <tjv, Tjv (<ts, 
Ts) are hermitian. We assume that the coupling matri- 
ces are diagonal in the quasiparticle e, h space: 



Wn,s = 



W] 



N,S 









(Al) 



where W^^^-j ~ ^(W^^(s))*- Since the Green's function 
of the (isolated) normal lead is also diagonal in the quasi- 
particle space, so will be Sjv — Diag(cr^ — iT%,(T 



• N 



iT}f), too. The self energy coming from the coupling to 
the superconductor has both diagonal and off-diagonal 
parts, but for £' < A it reads 



Sc = 



—he —h 



(A2) 



i.e. the superconducting lead does not broaden the lev- 
els. Moreover, since I/2) is only coupled with |/i) but 
not with any of the leads, the self-energy matrix ele- 
ments of the matrices Sat, will only affect the reso- 
nance energy £1 of the backbone state but not the energy 
£2 — £2 — Ep of the side coupled state. We now introduce 
the following notations: e^^ — £1 — Ep — {cr'^'^ + (J^g^) 
[where cf%^g are the (only) nonzero element of the ma- 
trices (T%% E^-" = {el'^ + e2)/2, Se'^'^ = (ef''' - £2)/2, 
= ± yf[5£e,h)2 j^^2^_ Denoting by Gbb{E) the 
retarded Green's function of the molecule and using the 
formulaM 



This work is supported partly by European Commis- 
sion Contract No. MRTN-CT-2003-504574 and by EP- 
SRC. 



RA^'Tr[nGBB{E)T%G\s{E)] 



(A3) 



to calculate the probability of the Andreev reflection, we 
find after straightforward calculations that 



APPENDIX A 

There are numerous equivalent approaches to calculate 
transport coefficients through phase coherent normal- 



Ra - 



AV%T%{E~£2f{E + e2f{a%^f 



Here the denominator is 



(A4) 
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D=[iE- e%){E - el) + i{E - e2)T%\[{E + eX){E + e^l) + i{E + e2)T%] - [affiE - e2){E + 82). 



(A5) 



and (tI'*, F^, are the only non-zero elements of the matrices cr|.'' ~ cr'i'^, F^, F^^. 
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